Analytical investigation of revival phenomena in periodically driven system 
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We present theoretical study of revival phenomena for a wave packet initially well localized in a 
one-dimensional potential in the presence of an external periodic modulating field. The classical 
motion, revival, and super-revival time scale are derived exactly for wave packet excited in the one 
dimension box potential. 
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The phenomena of collapse and revivals [1,2] of wave 
packet has been investigated both theoretically and ex- 
perimentally in a variety of physical system like that of 
Rydberg wave packets [3-5] , wave packets in semiconduc- 
tor quantum wells [6,7], coherent matter- wave in Bose- 
Einstein Condensation (BECs) [8, 9], and molecular sys- 
tems [10,11]. These collapse and revivals occurred due 
to the time evolution of wave packet driven by discrete 
eigenvalue spectrum and depend on the quantum number 
n. For the physical system with a non-quadratic depen- 
dence on quantum number n reveal a new sequence of 
collapse and revivals of the initially well localized wave 
packet. These are the super revivals [12] and is being 
observed in different physical systems [13]. 

In this article we present the time evolution and re- 
vivals of the wave packet in an arbitrary one-dimension 
potential in the presence of an external periodic field. 
The analytical analysis of this problem is explicitly dis- 
cuss in [14,15]. In these articles the authors discuss 
the problem by truncating the eigenvalue quadratically. 
Keeping the energy term up to second order in the Tay- 
lor series expansion, they have mapped their initially un- 
solvable time-depended Schrodinger equation into known 
Mathieu equation whose solution are Floquet solution. 
However, the mention approached is valid for small mod- 
ulation strength, i. e., A < 1. Moreover, the Mathieu 
characteristic parameter are expanded up to second or- 
der while ignoring the reaming term under the approxi- 
mation of small modulation strength. Here we consider 
the same problem by taking into account the energy term 
up to third order in the Taylor series expansion. Upon 
substituting into the Schrodinger equation, we obtained 
a third order differential equation. The solution of such 
equation is unknown. However, considering the term in 
the third derivative as a damping force, we again mapped 
our originally third order differential equation into known 
Mathieu-type equation. The procedure is valid and al- 
ready discussed in [16]. 

We denote the unperturbed energy eigenstate and 
eigenvalues of the correspondent time-independent sys- 
tem by \n) and E n , respectively, so that Ho\n) = E n \n), 



where Hq is the unperturbed Hamiltonian. The Hamil- 
tonian of the driven system in dimensionless form may 
be express as 



H = H + XV (x) sint 



(1) 



where A is the modulation strength and V(x) define the 
coupling. 

To obtained the solution for ^(i)) to the time- 
dependent Schrodinger equation we used the ansatz [14] , 
give by 

\m) = E C n (t)\n)exp{ i [E r + (n r)^]*-} (2) 

n 

where E r is the average energy of the wave packet in the 
iVth resonance. Substituting the above equation into the 
Schrodinger wave equation and project the sate vector 
(m\ onto the result, we have 

ihC m = E m - E r + (m - r)jj C m (t) + 777 y~V m ,n 

n 

x ^(n-m+N)t/N _ e -i{n-m+N)t/N j ^ ^ (3) 

where V m<n = (n\V(x)\m) are the matrix elements of 
coupling V(x). Under the approximation of keeping the 
stationary term m = n =fc N and dropping the fast os- 
cillating term, we obtained separate sets of coefficients. 
Moreover, assuming the matrix element independent of 
n, we can approximate by a constant V, i.e., V m ,m+N ~ 
V m , m -N = V. Hence we get N decoupled sets of equa- 
tion: 

ihC m = E m -E r +(m~r)^- Cm{t)+^-.{C m +N-C m ~N) 
L iv J 2,1 



(4) 



Expanding the energy up to n terms , we have 
ihC m = ( m -r)(E' r -^)C m + (^(m-r) 2 E' r 



+^(m - r) n E<^)C m + ^ Cm + N - Cm ~ N ) 

(5) 

we can write the above equation in angel representation 
by using the Fourier representation of C m as 

I-21T 



*Email: shahzadma@gmail.com 



c - - 



—i(r, 



2 



= 2N^l 9(O)e-^ m ' r)0/N d9 (6) 

with this particular choice of the Fourier representation, 
we obtained 



ihg(0,t) = Hg(0,t) 



(7) 



where 



H(0) = 



(i) n N n Er (>) d n 

(8) 



since the Hamiltonian H is time-independent, we can 
write the time evolution of g(0,t) as 



ff(M) = e(0)exp(p 



(9) 



using the above equation and considering the derivative 
up to third order, we have 

f iN 3 E'" rf 3 9(fl) N 3 E" d 2 e(9) , d®{6) 

{— — dB^ + — — &--< NE - h) ^r 

+ (-e + \Vsm(6))6{6)^ =0 

(10) 

Let 

?n 3 e'" n 3 e" , 

«=^f-, ? = 1 = i{NE'-h) {11) 

Eq. 10 can be rewritten as 



de 3 ' de 2 ' de 



-f—^ + (-e+\Vsm(9))e(8) = 

(12) 

Let us suppose the term in third derivative is the 
damping force. Then we can assume the driven force 
to be harmonic of angular frequency w[16]. Hence we 
can used the approximation 



de y 2(3 ) 



Or 



d 2 e 
de 



; 2 9 



(13) 



(14) 



in the damping term. So that Eq. 12 becomes 



- (- 2 + 7) ^ + ( - e + XV sin W ) 6 (0) = 

(15) 

Using the transformation 

e(^F(fl)exp(^(a^ 2 + 7 )0) (16) 



we have from Eq.(15) 

t{e)-^[( a uj 2 +1 ) 2 + Ap(e + \Vsin(e)))Y{e) = 



(17) 



by change of variable e = 2z + tt/2, we obtained a Math- 
ieu equation 

rW + {-4^(^ 2 +7) 2 -| + ^cos(2z)}F(z) = 



(18) 



this is Mathieu-type equation with the quantity 
— (l/4/3 2 )(a(x> 2 +7) 2 — {e/ f3) acts as the Mathieu constant 
a and the quantity (XV/ 0) acts as the Mathieu constant 
?■ 

Hence under the approximation of considering the 
third derivative as a damping force we have mapped our 
equation into know Mathieu equation. The Floquet so- 
lution of the Mathieu equation [17] may be written as 
Fj,(z) = exp(izv)P„(z) 7 where P v (z) are the Mathieu 
function [18]. 

Thus, the Floquet solution |\I> K (£)) to the Schrodinger 
equation can be written as 



\tp K (t)) =exp(ie K t/K)\u K (t)) 
where e K and \u K (t)) are given by 

£k = --^(auj 2 + -f) 2 - a K (3 
\u K (t)) = -^-^2exp(int/N) 



(19) 



(20) 



d<f>e^ N f^e^ n - r ^P v[k) \n) 

Jo 



(21) 



We calculate the autocorrelation function between the 
initial wave packet and the wave packet after certain evo- 
lution time t to get the super revival time. Since the 
wave function given in Eq.(19) form a complete set of 
basis vector, therefore we can write 



(m\*(t)) = J2\^\ 2et£ " t/h 



(22) 



where describe the probability amplitude in the nth 
state. Expanding the quasi-energy of the system around 
the resonant level r, the expression for the autocorrela- 
tion function becomes, 

(V(0)|*(t)> = ^|^| 2 ex P {-iM°) + (n-rV« 

n 

+ (n-r)V 2 > + -..]*} (23) 
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where lj^) denotes the jth derivative of e n with respect 
to n at n — r. From the above expression one can readily 
obtained the classical and quantum period for the driven 
system, i.e., 



Trl = 



2tt 



2-Kh 

on \n—r 



T A = 



2tt 



4-rrh 

o n \n—r 



(24) 



moreover, the super revival times appear in correspon- 
dence with the cubic terms of Taylor series expansion, 
that is 



2tt 

(3) 



CO 



12nh 

On \n— r 



(25) 



Considering the modulation term to be small, i.e., 
A <C 1, we can approximate the Mathieu characteristic 
parameter a v up to second order term in q. Therefore, 
Eq.(20) than becomes 



(26) 



taking into account the non-resonant situation, we have 
the explicitly relation for the classically revival time, that 
is 



T c i = T 



2v r 



-2v T 



2 (i/2 - 1)2 



(27) 



where T i = 2-Kh/ (3. Similarly we can obtained the ex- 
pression for quantum revival, 



T x = 2T 



2 + 



q 2 2(3^ r 2 - 1) 



2 - 1)3 



(28) 



Further, taking into account the cubic term of Taylor 
series we get super revival, 



T sr = 6Tq 



q Z 36^ 
2 (i/2 - l) 4 



(29) 



with this approached we analytically obtained the revival 
time in the presence of the an external periodically driven 
system. 

In summery, we examine explicitly the revival fea- 
ture in wave packet dynamics of a particle confined in a 
one-dimensional potential in the presence of an external 
periodic modulation modulated held. We have derived 
a rather complicated but accurate relation that reveals 
the super revival time as a function of the modulated 
strength. The calculation presented in this paper are 
more precise because of considering third power of quan- 
tum number dependence of the energy eigenvalue. 
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